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Nyquist  Frequency 


Abstract 


Nyquist  frequency  is  related  to  aliasing  and  the 
sampling  theorem.  This  report  explains  the  relation¬ 
ship  among  these  terms. 
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The  Nyquist  frequency  is  half  the  sampling  frequency  when 
a  continuous-tine  function  is  sampled  at  equally  spaced  tine 


points.  That  is,  the  Nyejuist  frequency  is (W&)  tin  radians 
per  unit  tine)  where  4  is^tjxe  tine  interval  between  two 

h'S  /  a nx>r*'  -  /  ^  e  u  S3*  s  * 

successive  sampled  data.  In-  this  article,  we— diecues  the 

ftp's."-*  It  <-  p 

basic  ideas  of  the  Nyquist  frequency  and  a  relevant  and  very 


P - important-theorem  «•  "the  sampling  theorem. 
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In  many  applications  involving  processing  a  continuous- time 
signal,  it  is  often  preferable  to  convert  the  continuous-time 
signal  to  a  discrete-time  signal  since  discrete- time  signal 
processing  can  be  implemented  with  a  digital  computer.  It  is 
important  to  examine  whether  the  discrete-time  signal  preserves 
all  the  information  in  the  original  continuous- time  signal. 

We  first  consider  the  case  that  the  signal  x(t)  is  a  real-valued 
function.  Assume  that  its  Fourier  transform  X(<d) 

CO  oo 

=  J*  x(t)exp(-iut)dt  exists  and  that  x(t)  =  ^  f  X(u))exp(iwt)  dc 

—CO  —CO 

for  all  t.  The  signal  x(t)  is  sampled  at  t  =  nA,  n  =  ..., 
-1,0,1,...  .  We  are  interested  in  interpolating  x(t)  from 
its  samples  x(nA).  A  natural  question  arises:  Under  what 
conditions  can  x(t)  be  perfectly  reconstructed  from  x(nA)? 

The  samples  x(nA)  can  be  related  to  X(w)  as  follows. 


x(nA)  =  yjf  S  XU>)exp(iwnA)dto 


•>  ir/A 

=  s-  f  m  X(w  +  2irk/A)  ] exp CiwnA) dw 
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Therefore ,  X^Cco)  =  ZX(w  +  2irk/A)  C|  w|  <  w/A)  is  the  discrete 
Fourier  transform  of  the  sequence  x(nA).  Obviously,  xd(*)  is 
obtained  by  folding  X ( * )  every  tt/A  radians  per  unit  time. 

(.Here  we  identify  with  -w.)  This  frequency  tt/A  is  the 
Nyquist  frequency,  and  hence  also  called  the  folding  frequency. 
It  is  easy  to  see  that  X(*)  is  not  uniquely  determined  by  Xd(*). 
In  other  words,  some  sinusoidal  components  of  different 
frequencies  (e.g.  2TTk/A+w0,  k  *  ...,-1,0,1,...)  in  x(t)  cannot 
be  distinguished  from  one  another  by  the  observations  x(nA) . 
This  is  called  aliasing.  Aliasing  is  the  effect  of  under¬ 
sampling.  This  effect  is  the  principle  on  which  the  strobo¬ 
scopic  effect  is  based  [6,  Section  8.3]. 

When  x(t)  is  a  band-limited  signal  with  X(w)  *  0  for 
|  h)|  >  <*>M/Xd(u>)  is  identical  to  X(w)  (i.e.  no  aliasing)  if 
ojm  <  tt/A.  In  other  words,  from  the  uniqueness  property  of 
Fourier  transform,  x(t)  is  uniquely  determined  by  its  samples 
x(nA)  under  the  condition  that  X(u>)  =  0  for  |u>|>  tt/A.  This 
is  usually  called  the  (Shannon)  sampling  theorem  on  infor¬ 
mation  theory  [5].  From  the  sampling  theorem,  if  we  sample 
the  signal  x(t)  at  a  rate  at  least  twice  the  highest  frequency 
in  x(t),  then  x(t)  can  be  completely  recovered  from  the 
samples.  This  sampling  rate  (twice  the  highest  frequency  in 
x(t))  is  commonly  referred  to  as  the  Nyquist  rate.  Actually 
x(t)  can  be  explicitly  written,  in  terms  of  x(nA) ,  as 


It  should  be  noted  that  band-limited  signals  are  generally  not 

realizable  physically,  for  xlt)  *  ^  f  X  (.<»>)  exp  (iwt)du  is 

*  ^ 

analytic  in  t,  as  a  complex  variable,  and  therefore  cannot 
vanish  for  all  t<-T  for  arbitrarily  large  T.  Therefore, 
aliasing  is  inevitable  in  practice.  A  discussion  on  error 
bounds  for  aliasing  can  be  found  in  15] . 

In  some  applications,  the  signal  xtt)  is  assumed  to  be 
bandpass,  i.e.  there  exist  0  _<  <  uj1  such  that  X(w)  =  0 

outside  the  intervals  Iwq,w11  and  [-w1,-o)q].  The  sampling 
theorem  says  that  x(t)  can  be  recovered  from  equally  spaced 
sampling  at  a  rate  of  2uj^.  Actually,  this  rate  20^  is  too 
conservative.  It  has  been  shown  [4,  Section  8.5]  that  a 
sampling  rate  of  is  enough  to  recover  x(t)  where  v 

is  the  largest  integer  not  beyond  (u^-uig)  . 

The  sampling  theorem  has  been  generalized  to  many  situations 
such  as  random  signals.  When  x(t)  (-«<t<°°)  is  a  wide-sense 
stationary  stochastic  process,  possessing  a  spectral  density 
which  vanishes  outside  the  interval  [— tt/A , it /A],  Balakrishnan 
showed  [1]  that  x(t)  has  the  representation 


x(t) 


N 

lim  Z  xCnA) 
N  »  n=  -N 


sin  it  { Ct/A)  -n) 
it  { Ct/A) -n} 


for  every  t,  where  lim  stands  for  limit  in  the  mean  square. 
Gardner  [3]  derived  a  similar  result  for  non-s tat ionary 
stochastic  processes.  Obviously,  the  Nyquist  frequency  and 
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Nyguist  rate  can  be  similarly  defined  in.  the  random  signal 
case. 

Blackman  and  Tukey  12,  Section  12]  provided  good  inter- 
pretations  on  aliasing.  Jerri  I 5 J  gave  an  excellent  review 
of  the  sampling  theorem  and  its  various  extensions  and 
applications.  He  discussed  topics  such  as  unequally  spaced 
sampling,  higher-dimensional  functions,  nonband-limited 
functions  and  error  bounds  for  the  truncation,  aliasing, 
and  jitter.  A  very  exhaustive  bibliography  can  be  found 
therein. 
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